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^^«. ' Abstract. We consider an inverse scattering problem for Schrodinger operators 

Q>>. ' with energy dependent potentials. The inverse problem is formulated as a Riemann- 

Qs , Hilbert problem on a Riemann surface. A vanishing lemma is proved for two distinct 

symmetry classes. As an application we prove global existence theorems for the two 

distinct systems of partial differential equations ut + [y? 1'^ + w)^; =0, uit ± Uxxx + 

>^ ' (uw)x = for suitably restricted, complementary classes of initial data. 

<D ■ 

Q 
o\ ■ 

Dedicated to Hugh Turrittin on his 90th birthday 

1. Introduction 
In this paper the scattering theory of the energy dependent Schrodinger operator 

iD^ + P)ij^{ikp{x)+q{x))^P, D = -^, (1.1) 

ax 

0^ , is investigated and used to prove global existence theorems for the "isospectral 

flows" : 
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Mt + (w; + — )^ = 0, Wt + Urcxx + {wu)x ^0; 
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^' Ut+(w+—)a;^0, Wt-Uxxx+{wu)x ^0 
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The scattering problem for such energy dependent Schrodinger operators was 
first considered by Jaulent [8], Jaulent and Jean [9], Kaup [10], and more recently 
by the present authors [14]. In [8], [9] and [14] the problem was considered for 
purely imaginary p and real q tending to zero at infinity. The inverse problem was 
solved by a modification of the Gel'fand-Levitan-Marchenko method. 

Kaup studied the scattering problem under the assumption of real p, and q 
tending to non-zero limits at infinity, in connection with a long wave approximation 
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of Boussincsq type (equations (1.2)). He obtained a coupled pair of equations of 
GLM type, but did not investigate their solvability. 

The GLM equations are of Fredholm type, and so their solvability is a conse- 
quence of uniqueness. Such uniqueness theorems are valid, (though usually not 
proved in the literature) in the case of the standard Schrodinger equation, or the 
isospcctral problem obtained by Zakharov and Shabat in their study of the non- 
linear Schrodinger equation [15]. But they fail in general for 2 x 2 AKNS systems 
without further symmetry assumptions on the potential, for example, that it is 
symmetric or skew symmetric, Hermitian or skew-Hermitian. 

We investigate the solvability of the inverse scattering problem for (1.1) in the two 
cases p real or imaginary, corresponding to the two flows (1.2) and (1.3). Uniqueness 
theorems are proved in these two cases under complementary, restricted conditions 
on the scattering data. Global existence theorems for the flows (1.2), (1.3) can 
then be proved by showing that the constraints on the scattering data are invariant 
under the flows. 

This, however, entails restrictions on the initial data. Since global existence 
theorems can be proved for the two flows only for restricted initial data, the phys- 
ical relevance of these two models is unclear. Our interest in the matter stems 
rather from the novel features of the associated inverse scattering problem, and its 
formulation as a Riemann-Hilbert problem on a Riemann surface. 

Riemann-Hilbert problems arise in integrable systems as a more general formula- 
tion of inverse scattering problems. In second order cases, such as the Schrodinger 
and Zakharov-Shabat scattering problems, the Gel'fand-Levitan-Marchenko theory 
and Riemann-Hilbert method are equivalent, being essentially Fourier transforms of 
one another. But there is no GLM theory for nx n first order systems or n order 
ordinary differential operators (n > 2); and the inverse scattering problem must be 
formulated as a Riemann-Hilbert problem [1, 2]. Moreover, even for second order 
problems the method of Riemann-Hilbert problems, coupled with the method of 
steepest descent, can be used to obtain precise asymptotic behavior of the solutions 
[5]. 

One can go from (1.2) to (1.3) by the transformation 

u{x, t) -^ iu{x, it), w{x, t) — > —w{x, it)] 

but this transformation is complex, so the two equations may be considered to be 
distinct real forms of one another. 

Sachs [13] found a pair of tau functions for (1.2) and used them to construct 
rational solutions similar to the Calogero-Moser solutions of the KdV equation. 
Matveev and Yavor [11] used Riemann surface theory to construct finite gap solu- 
tions, and, in the limiting case, multi-soliton solutions. 



Isospectral flows for energy dependent operators arise in a wide variety of ap- 
plications. Recently, for example, Camassa and Holm [4] found a special shallow 
water wave equation generated by an energy dependent iso-spectral operator re- 
lated to the Dym hierarchy. For a formal discussion of isospectral flows generated 
by operators with energy dependent potentials, see Fordy [7]. 

The Gel'fand-Dikki hierarchy of isospectral flows of n*'* order ordinary differen- 
tial operators can be formally extended to operators L in which the coefhcient of 
D^ is a polynomial of degree j — 1 in the spectral parameter. These considerations 
suggest, at least on a formal level, the existence of a large class of flows, of which the 
Gel'fand-Dikki flows are a special case. We conjecture that all these flows are local. 
The corresponding inverse problems lead to a broader class of Riemann-Hilbert 
problems than those considered in [3]. 

2. The forward problem 
Throughout this paper we assume that q ^t 1 as x ^t ±oo and rewrite (1.1) as 

{D^ + E^ -{ikp + q))i, = Q, E^ = k^ + l, (2.1) 

where now p and q belong to Schwartz class S. Although the restriction to potentials 
in Schwartz class could be weakened considerably in the theory of the forward and 
inverse scattering problem, the class is a convenient one for discussing completely 
integrable systems with an infinite number of conservation laws, and so we shall 
make this assumption throughout. 

The scattering problem may be formulated as a Riemann-Hilbert problem on 
the Riemann surface E^ = k'^ + 1. Following Kaup, we introduce the uniformizing 
parameter z: 

The transformation z i-^ E is a double covering of the E plane by the z plane. 
Define four regions Uj, j = 1, 2, 3, 4, where ZYi = {z : 3 z > n |z| > 1}, ZY2 is its 
image under complex conjugation, U3 the upper half of the unit disk, and Ui the 
lower half of the unit disk. Let 

U+=UiUU4^{z:^E> 0}, U-=U2UU3 = {z:^E < 0}. 

Finally, let S denote the union of the real line and the unit circle in the z plane, 
less the origin. 

We construct two sets of wave functions oi (2.1), tp and (p defined by the following 
asymptotic behavior: 

V'±(x,z)=m5;^*e±'-^^, zeU± </)±(x, z) = m!ge^'-^^, zeU± 



where to™* and to,™ are analytic on IA± respectively, and 



lini ?Ti'^ =: lira to,^" = 1. 

X — ^ — oo X — *oo 

Equation (2.1) can be converted to a Volterra integral equation, for example, 
V'±(x,z) = e±'^- ^1°° ^J^I^^zA^^kp{y) + q{y))^±{y,z)dy, zeU+. 
This leads to the following integral equations for to.™*: 

mT{x,z)^lTJ 2lB {^kp{y)+q{y))m°^\y,z)dy, zeU±. 

(2.3a) 
The wave functions (j)± and m™ satisfy similar Volterra integral equations. For 
example, 

w}1{x,z) = ItI ^~%.^ -{ikp{y) + q{y))m'^{y,z)dy, z€ZY±. (2.3b) 

J-oo 2iii 

The integral equations for m™, to-™* may be solved by successive approximations 
for all z e U± \ {0, oo}. (See also the form of the wave functions given in the 
Appendix.) 

Throughout this paper we shall understand by ^ the scctionally analytic function 

rV+(a;,z) z£U+; 
\ ip-{x, z) z £ U-. 

The other wave functions are denoted similarly. 

The asymptotic behaviors of the wave functions to.™*, m}^ as 2; — > 0, oo are 
easily determined. We note that 

k f 1 z ^ 00; 

To simplify the discussion, we assume, as in [14], that 

pdy — Q. 



This assumption is not essential and could be dropped, though with a consequent 
increase in the computational details. Letting 

f-OC 

\ = ePI\ P= / p{y)dy, 



we easily deduce from the integral equations (2.3a,b), that 

lim mrVx,z)=AT\ lim m5y*(x, z) = A±\ z^U±] 

z— ►oo z^O 

(2.4) 
lim m'?.(x, z) =A^"^, lim m'P(x, z) = X^^ z eU±. 

z~*oc z^O 

Note that A is real when p is real and has modulus 1 when p is imaginary. 

If p and q are in the Schwartz class, the reduced wave functions m°"* have 
asymptotic expansions in powers of z. Substituting tp+ = m°^*^e'^^^ into (2.1) we 
get the following equation for m — tti™*: 

m + 2iEm — (ikp + q)m = 0, 

We seek a formal asymptotic expansion of m in inverse powers of z; that is, we set 

oo 

m{x,z)^y mj{x)z^-' 
j=o 

and substitute this asymptotic expansion into the differential equation for m. At 
orders and 1 we obtain 

zttIq — -pniQ — 0, ttIq + inii qniQ = U. 

iFrora these two equations we obtain 

mo = e-^/2 = X-\ , = ^ + ^ (I + ^m,e//^) (2.5) 

As in the theory of the standard Schrodingcr operator, the wave functions satisfy 
the usual relations on S: 

(2.6) 

cj,^{x,0=c{Oi^-{x,0+d{Oi^+{x,0 

where (/)+(x,^) denotes the limiting value of (j)^{x,z) as z ^ ^ from ZY+, etc. The 
coefficients in (2.6) relate the incoming and scattered waves. Thus, for example, 
from (2.6) it follows that for ^ G S 

r e^*-^^ X -^ — oo; 

'^+^'''^^^ la(C)e-*^-+6(C)e*^^ x ^ oo. 

The matrix of coefficients 

\ c a 



is sometimes called the scattering matrix ([6]). [?].) The coefBcients of S may be 
computed in terms of Wronskians of the wave functions, just as for the standard 
Schrodinger equation. Since the Wronskian of solutions of (2.1) is independent of 
X, we have VF (</)_, (/)+) = VF(e'-^^, e"*-^"=) = 2iE. It then follows from (2.6) that, 

"(^) ^ 2^E ' ^(^) = —^2l^^ 

(2.7) 

'^^^ ^ 2^E ' '^(^^ = -2»g ■ 

Moreover, 

W{(l)+,(j)^) =-2iE = VK(aV'- + b^p+, ctp- + #+) = W{ip-,ij+){ad - 5c) 
==-2i£;det 5; 

hence det S = I. 

At a zero of a{z) at z^ G Z^+ the wave functions </)+ and ^+ are linearly depen- 
dent, and we can write 



9± 



{x,zj) ^cjip±{x,z+). 



where c^ are the associated coupling coefficients. Since (f>^ decays exponentially as 
a; ^ — oo and ?/;+ decays exponentially as a; ^ oo, this common function belongs 
to L2 hence constitutes an eigenfunction of (2.1), called the bound state. Similar 
considerations hold on W_. Thus, the bound states of the problem are given by 
zeroes of a{z) in U^ and d{z) in U^. We denote these zeroes by z respectively, 
and the associated coupling coefficients by c ■ . 
We define generalized reflection coefficients 

When the zeroes of a and d are simple, the scattering data for the operator (2.1) is 
the set 

{r±{0, eeS; zf ; cf } 

Under the mapping z h^ — 1/z, E h^ —E and fc h^ fc so the differential equation 
(2.1) is invariant. We denote this mapping by s, and write s.tp{x, z) = ^{x, — 1/z). 

Lemma 2.1. Under the mapping s the wave functions transform as s.ijj+ = ip-, 
s.jTi^"' = TO°"*, i.e. '0+(x,— 1/z) = ^_(x,z), etc. Moreover, for z G U^, a{z) = 
d{-l/z); while h{£,) = c(-l/^) and r+{i) = r_(-l/0 for ^ (^ E. 

Proof. Under the mapping z h^ ^1/^ the Voltcrra integral equations in (2.3a) for 
^* are interchanged. Since the solutions are uniquely determined, they 



are the same. The result extends immediately to ')p± and (j)± . The relations on the 
generalized reflection coefficients then follow from their expressions in terms of the 
Wronskians. D 

When q is real and p is real or imaginary, the wave functions possess additional 



symmetries. Since k — \/E^ — 1, k is real on E^ > 1 and imaginary on the slit 
— 1 < i? < 1. Hence (2.1) is invariant under Schwarz reflection across the slit Up 
is real, and under Schwarz reflection across E'^ > 1 if p is imaginary. The slit lifts 
to the unit circle in the z-plane, while the rays E'^ > 1 lift to the real line in the z 
plane. These observations lead to the following result: 

Lemma 2.2. When p is real, (2.1) is invariant under Schwarz reflection across 
the unit circle and across the imaginary axis in the z -plane; and the wave functions 
(j) and -0 satisfy the symmetry conditions (/'+(x, z) = 4>-{x, 1/z) = </'+(a;, —z). 

When p is imaginary, (2. 1) is invariant under Schwarz reflection across the real 
axis 3 z = 0, as well as under the reflection z -^ —z^^ ; the wave functions possess 
the corresponding symmetries: (l)+{x,z) — 0_(x, f) = 0+(x, — 1/z). 

Proof. Take the case of real p, and consider the wave function (f>. An easy calculation 
shows that, for z <E Z^+, both (t>-{x, 1/z), and (p+ix, z) are asymptotic to e~*^^ as 
X -^ — oo and satisfy the same differential equation. Since solutions of (2.1) are 
uniquely determined by their asymptotics as x ^ — oo, they in fact coincide. The 
other cases are handled in the same way. D 

Theorem 2.3. When p is real or imaginary the scattering data possess additional 
symmetries: 

For p real, a(^) — a(— C) — "^(l/O '^'^'^ ^+(C) ~ i~+{^0 = ''-(I/O- Moreover, if 
Zj is a bound state, then so are —'zj, 'zj^^ and —z^ . The four associated coupling 
coefficients are respectively Cj , cj, cj, and Cj . 



When p is imaginary, a(^) = a(— 1/^) = rf(^) and r+(^) = r+(— 1/^) = r_(^). 
The bound states also appear in fours: Zj, zJ, — z7 , — "zj^^, with corresponding 
coupling coefficients Cj, cJ, Cj, cJ. 

Finally, 

I ^4 = L) p imaginary. 

Proof. The lemma follows from Lemma 2.2 and the computation of the coefficients 
of S in terms of Wronskians. /^From det S ^ ad — be ^ 1 we get 1 — r-^-r- = 1/ ad] 
and the last statement above follows from the relationships between a and d, r_|_ 
and r_ on S when p is real or imaginary. D 

The following result will be useful in our discussion of the vanishing lemma, in 



Lemma 2.4. Let p he real and consider the hound states lying on iR \ {0}. Then 

1 f°° , .of iuJia'(iu!j)cj 1 < w," < oo; 

-/ (2 cosh a, ~p)<l>',dx^\ '\;.l n^\.i (^Ja) 

^ J_oo I —iOJjd (lOJjjCj {) < UJj < 1, 

and 

1 f^ \ 2 \ —iujja'{iujj)cj —1 < LOj < 0; 



(2cosha,+p)0^^dx= , ': '" /' 1 (2.9b) 

^ J_oo I lUJjd [lUJjjCj — oo < tJj < —1, 

where \ujj\ ~ e"^ m 6ot/i cases, Cj is the coupling coefficient associated with the 
hound state, and (j}j ~ (f>±{x,iLUj) according as luj G U±. 

Proof. Differentiating (2.1) with respect to z we have 

(1)2 + E^ - V)(t>' = [V - 2EE')(t> (L»2 + ^2 - V)(t) = 0, 

where primes denote differentiation with respect to z, and V = ikp + q. Multiplying 
the first of these equations by 0, the second by (j)' and subtracting, we obtain 

DW{(j), (f)') = [V - 2EE')(j}'^, W{(j), (j)') = (t)D(j)' - (P'Dcf). 

Now (j)j — 4>+{x, iuj) tends to zero exponentially as x — > ±cx), while (/>'■ tends to zero 
exponentially as x -^ — oo. Integrating the above expression over the real line we 
obtain 

i f°° 
lim W{(f>j, (/)' ) — — / [piy) sinh Uj — 2 sinh a,- cosh aM^ dy 

/OO 
[piy) -2 cosh aj](j>'^jdy, 
-oo 



i sinh a ''°° 



where loj = e"^ for the case luj > 0. ^From (2.7), W{4'+,'ip+) = 2iEa{z) for 
z G W+; differentiating this identity with respect to z, and evaluating at a bound 
state z = iujj we obtain 

W{(j)'j,ipj) + W{(j).j,ip'j) = 2iEja'{iujj) = -2 sinh aja'{iujj), 

with (j)': = (j)'^{x,iujj), etc. 

We now use the relation (/)j = c&i'j: in the above result, and let a; ^ oo. The 
left side becomes 

The second term tends to zero, since V'j ^nd i/^' tend to zero exponentially as 
X -^ oo. The identity (2.9a) now follows. The identity (2.9b) is obtained in the 
same manner. D 

As an immediate consequence of Lemma 2.4 we have: 



Theorem 2.5. Let p be real and \p\ < 2. Then all the bound states Zj — iujj G 
iK. \ {0} are simple, and Cj > 0, where 



--^ 1 < ujj < oo; 



c. 



Ci = < a"' " A] = 17^ (2-10) 

We now formulate the forward problem as a Ricmann-Hilbcrt problem. Define 
the row vector m by 



(2.11) 



?n(x, z) = i 

Here, as usual, CT3 denotes the Pauli spin matrix 

A simple computation shows that (2.6) are equivalent to the jump relations 



1 -r+e2*^^ \ , ^.( 1 

1 =™-(^'0 _,_e-2.B- 1 



We write this as 

m+{x,0=ra-{x,i)v{x,0, CeS (2.12) 

where 

<-'^)-(-._(Je— Yi^l'Z)- ^'-''^ 

The row vector tti has prescribed asymptotic behavior as a; ^ 00 and as z tends 
to or 00, namely: 

?n(a;, z) ^ (1, 1), as a; ^ 00 (2-14) 

and by (2.4) 

r (A,A-i) z^O; 
m(x,z)^ '1 , (2.15) 

1^ (A , A) z ^ 00. 

Moreover, 

^+(a;,i) = m'["*(a;,i) =-(/>_ (a;,i) =m°^\x,i). (2.16) 

This identity follows by observing that i/'±(a;,i) satisfy the same Volterra integral 
equation, viz. 

^±{x,i) = 1- {x-y){q- p)i}}± (y, i) dy. 

J X 

We also have iJ-{x, —i) = tp+ix, —i). 

In the case of a reflcctionless potential, r^{S,) = 0, and vix^S,) — I everywhere 
on S. 
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Theorem 2.6. The inequality 1^(^)1^ < 1 holds on the entire real line in the case p 
imaginary, and everywhere on the unit circle except at ±i in the case p real. When 
p is real we have generically, (that is, except in the case of reflectionless potentials) 
r+(±i) = r_(±i) = — 1 and 

v{x,±i) = L ^ 

In all cases, r±(0) = 0, and v{x,0) = I. 

Proof. By Theorem 2.3 we have |ap — |6p = 1 on cither the real hne or the unit cir- 
cle, according as p is imaginary or real; hence a cannot vanish on the corresponding 
portion of E. Writing this identity as 1 — |t'('C)P = k(C)l~^7 "^^ see that |r(^)| < 1 

whenever 

1 _ 2iE 

^} ^ w {<!>+, ij+) ^ • 

We see that this holds whenever E y^ 0, hence everywhere except at ^ = ±i. 
Writing the first equation in (2.6) in the form 

^/^'x == ^-ix,i) +r+{£,)ijj+{x,£,), 

and noting that l/a{i) — (note that, since we have already observed that a{i) ^ 
for real p, W{(j)^{x,i),'ip^{x,i)) ^ 0) we see that r^(i) = 0. A similar argument 
shows that r_(±«) = 0. The expression for w(x,±i) above follows from the fact 
that E{±i) = 0. 

We write r+(0 = -W^(<?!>+, V'-)/VK(0+, V'+)- Then 

As ^ ^ 0, toT* -> A, and mV" ^ A"\ so that 



d 



lim^VK(m![?,m^"*) ^W{\'^,\) = -2— log A = p(x) 



A similar calculation shows that 

lim W((/)+,V-) = Ymie-'^^^''W(wtl,m°^^) ==0 

since m™ and to'" both tend to A^^ as ^ ^ 0. Noting that E{^ ^ oo as ^ ^ 0, 
we see that r+(^) ^ as ^ ^ 0. A similar argument shows that r_(0) = 0. D 

Since f /a(«) ~ ^/d{i) = we have, in the generic case, 

TO+(x,i) = (TO7*(x,i), 0), TO_(a;,i) = (0, TO™'(x,i)); (2.17) 
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Lemma 2.7. The row vector m possesses the symmetries 

s.m^m{x,-l/z)^m{.,z)R, R ^ (\ J 



and v{x,^) = Rv ^{x,—llS,)R. When q is real, m{x,z) — m(x,l/ z)R if p is real; 
and m(x, z) — m{x, z)R when p is imaginary. 

The row vector m has poles at the zeroes of a and d. At a simple pole of m, 
there is a triangular matrix Vj{z) such that 

Pj{x,z) = me''^'"'^Vj{z)er''^'"'^, 

is regular at z = Zj [1]. In fact, from the bound state relation 4>±{x,z, ) = 
cfip±{x, z^) wc find this to be true if we take 

^'^From the symmetry a{z) = d{—l/z), it follows that if m has a pole at Zj G U+, 
with coupling coefficient Cj, then it also has one at — 1/zj G U-, with the same 
coupling coefficient. In fact, from (t>+{x, Zj) ~ CjiI)j^{x, Zj) and the s-symmetry we 
find (/)_(a;, — 1/zj) ~ CjXp^{x,—l/ Zj). ^From the identity a'(z) = — z^^(i'(— 1/z), 
we find 

Rl ^ £z ^ 4z. (2 19) 

(^l/^^j) ZjU'iZj) Zj 

^From (2.18) we deduce that mi is regular in U+ while TO2 is regular in ZY_, and 
that 

, , r (0,-A,TOi(x,z,)e2»=^^0, z,eU+; 

^^'"'^'^ == n / \ 2..iJ. ^^ ;, 2.20 

z=z, I {-Djm2ix,Zj)e ^*=^-^^0), ZjeU^, 

where 2Ej = (zj + l/zj). 

It is easily seen that both components of the row vector pe^^^'^^ satisfy the 

differential equation (2.1). Moreover, pj satisfies (2.14), (2.15). Now choose a small 

circle Cj around Zj, containing no other bound states and not intersecting E, and 

define a new vector m by 

Pj{x, z) z inside Cj] 



m{x,z) 

m{x, z) z outside Cj. 

Then the jump of m across the contour Cj is 

m+ = m_Wj(x,z), Vj{x,z) = e''^''"'''vj{z)e~''^''""- , (2.21) 

where the limits m± are taken relative to a counterclockwise orientation of Cj and 
Ej = E{zj). If there are only a finite number of simple poles, we can modify 
the original contour by adding a small circle around each pole, and redefining m 
accordingly, so that the original Riemann-Hilbert problem is replaced by a modified 
one on the extended contour. 
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Theorem 2.8. Let p, q (z S, and let m he the reduced vector valued wave function 
defined in (2.11). Then: 

(1) m satisfies the Riemann-Hilbert problem (2.12, 2.21) with v and Vj given 
in (2.13), (2.18); 

(2) m satisfies the asymptotic conditions (2.14), (2.15), and the conditions 
(2.16) at z — i together with corresponding conditions at z = —i; 

(3) in the generic case, m satisfies the constraints (2.17), plus similar con- 
straints at ^ ^ —i; 

(4) m satisfies the .symmetries of Lemma 2. 7. 

3. Two Vanishing Lemmas 

A general approach to the solution of Riemann-Hilbert problems such as (2.8) 
is given in [3] . The inverse problem is reduced to the inversion of a linear operator 
of the form L + Ti + T2 where Ti is of small norm and T2 is compact. Once 
this reduction has been obtained, the Fredholm alternative applies, and the proof 
of existence and uniqueness is reduced to proving a uniqueness theorem for the 
Riemann-Hilbert problem. Such a uniqueness theorem in the context of Riemann- 
Hilbert problems is sometimes called a "Vanishing Lemma" [3]. We prove the 
vanishing lemma for our problem in this section. The remaining details of the 
inverse problem arc the same as, but somewhat simpler than, those given in Part 
H of [3]. The inverse problem treated in that book, the Riemann-Hilbert problem 
for an n*'* order ordinary differential operator, is complicated by the intersection 
of rays at the origin. 

The proof of the vanishing lemma is based on the symmetries on the reflection 
coefficients given in Theorem 2.3. For the purposes of the vanishing lemma, x 
appears only as a parameter, so we suppress it for the discussion. We first prove a 
vanishing lemma for the Riemann-Hilbert problem that arises in the case when p 
is imaginary and there are no bound states. 

Theorem 3.1. Let m be a piecewise analytic function of z m 3z 7^ 0, with bound- 
ary values m± G L2(R); and let m satisfy the Riemann-Hilbert problem 

m_|_(^) = m_(^)u(^), ^ G R; m — 0{z^ ) as z ^ 00, 

where 



-(0=1 ^ . ']tLa]- (3.1) 



Assume \r\ < 1 a.e. Then m vanishes identically 



Proof. We define 
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where 

mJ — I 

1712 



Thus, m{z) ■ m^{z) — mi{z)mi{z) + m2{z)m2{z). 

By the Schwarz reflection principle, F is sectionally analytic in 3z ^ 0. Its 
boundary values on the real line are 

F+{0 ^ ^ ; 

and 

^-^^> 2 - 2 • 

Therefore, the jump of F across the real line is 

[F] =F+-F^=m^ (^^) "^- = l^-(^)l' + (1 - HOn^-m'- 

Since F has no poles and is 0(z^^) as z ^ oo, we can integrate it around a circle 
of radius _/?„ centered at the origin, and let _/?„ ^ oo to obtain 

[F] dC ^ / \ml{0\' + (1 - \riOr}\mUO\' d^ = 0. 



Since \r\ < 1 a.e. the boundary values m^{^) and m^{S^) = m^(^)v{^) vanish 
a.e. on the real line. Since to — > as z ^ oo and has no other singularities, to. 
must vanish identically. D 

We now prove a vanishing lemma for the case corresponding to real p. The 
operator (2.1) is self-adjoint only for z on the imaginary axis or the unit circle, so 
in this case we restrict the scattering data to lie on the union of the imaginary axis 
and the unit circle. Note also that the data and wave functions are invariant under 
Schwarz reflection in the unit circle and in the imaginary axis. 

Theorem 3.2. Let m{z) be a piecewise meromorphic function of z in \z\ ^ 1 
satisfying the Riemann-Hilbert problem to+(0 = m^{^)v{^) on \^\ = 1, where v 
has the form (3.1) on the unit circle. Assume that \r{S,)\ < 1 a.e. on the unit circle, 
and that 

(1) the poles of m are simple and lie on the imaginary axis; 

(2) the residue matrices at Zj — icuj and z^j — ij^j are given by (2.18a) and 
(2.18b) respectively, where Aj and Dj satisfy (2.19) and Aj/zj < 0. 

(3) both components of m are 0{z^^) as z ^ oo; 
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Then m vanishes identically. 

Proof. We proceed in three steps. 

Step 1. We first prove the theorem in the case there are no bound states. Define 
the function 

F(z) = { 2. I I , 

By the Schwarz reflection principle, F is sectionally analytic in \z\ ^ 1; and also 

F{z) = <^ , 3.2 

^ ^ 1 0(1) z->0. ' 

We orient the unit circle in the counterclockwise direction, and denote the lim- 
iting values of F from the interior and exterior of the unit circle by i^^ , _F_ respec- 
tively. On 1^1 = 1, 

F+{i) - ^^ - ^^ ; 

and 

p ,.. ^ m+(i) -toUO ^ m_(^)i;(g)m^g) 
" 2^ ^ 2^ ■ 

(Note that ^ = 1/^ on the unit circle.) The jump of F across the unit circle is 
therefore 

( V + V \ + 
[F\ =m^ I — ^T — I iTi 



2e 



lei = 1. (3.3) 



Let r^ be the contour consisting of the unit circle traversed in the clockwise 
direction and the circle \z\ ^ Ra > 1 in the counterclockwise direction; and let FJj 
consist of the unit circle traversed in the counterclockwise direction together with 
the circle \z\ ~ 1/Ra traversed in the clockwise direction. We let Tj^^ ~ Fjj UF^ . 
For all i?a > 1 we get, by Cauchy's theorem, 

since F has no poles. Letting Ra ^ oo we obtain, in virtue of (3.2), 

Now d^/i^ is a positive measure on the unit circle oriented in the counterclockwise 
direction, and |r_|_(^)| < 1 a.e. The boundary values ?Ti_ and m^ = m^v must 
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therefore vanish a.e. on the unit eircle. Since m has no singularities in the finite 
plane and vanishes at z — oo, m vanishes identically. This completes step 1. 

Step 2. We next prove the theorem when m has simple poles at iuj and i/coj, 
where ujj > 1 and 1 < j < A^ and there are no jumps across the unit circle. We 
define 



F{z) = 



M-y^H--^) \z\>i- 



2z 



|z|<l. 

We now calculate the residue of the first term at a pole Zj — iujj, ujj > 1. Since 
LOj > 1, mi is regular at Zj — iujj, and, by (2.20), 



Res?7i2(z) 



-Ajmi{iujj). 



Now observe that if the function f{z) has a simple pole at zq with residue /o, then 
the function ,f{—z) has a pole at — zo with residue — /q. Therefore 



Resm2(— z) 



Ajmi{iujj), 



and 



Res 



mi{z)m2{~ z) 



2z 



2iLU, 



\mi{iujj)\^ 



2z, 



\mi{iujj)\'' 



(Note that, since Aj/zj < 0, Aj is imaginary.) 

The residue of the second term is calculated in a similar manner, and we find 
that 



ResF 



A- 



A similar calculation yields 
ResF 



D 



A- 



z—i/iUj — f/^ 

Here we have used (2.19). 

We again integrate F around the contour Tji^ and let Ra ^ oo. By the Cauchy 
residue theorem, 

hm -i- / F{z)dz = 2j2C,\mlizu;,)\' = 0. 

where Cj = —Aj/zj. 

Since Cj > by assumption 2, each term in the above sum must vanish. Since m 
has no singularities in the plane and tends to at infinity, m must vanish identically. 
This completes step 2. 

Step 3. When in has both poles and jumps across the unit circle, we construct a 
gauge transformation that removes the poles and reduces the problem to the case 
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where there are only jumps. We construct a matrix valued function of z with the 
following properties: 

i) M is meromorphic in the extended complex plane with simple poles at Zj and 
M ^ yl as z ^ oo, where A is an invertible matrix. 

ii) Mvj is regular at Zj. 

iii) M(^) is unitary on the unit circle. 

Given such an M we put m ~ wM , where m is the null solution of the original 
problem. Since mvj and Mvj are regular at Zj, it follows that w is regular at Zj. 
Moreover, w = 0{z^^) at infinity. The jumps of w on the unit circle are given by 

w+ = w-MvM^^ = w-MvM* 

since M is unitary on the unit circle. We now apply the argument of step 1 to the 
row vector w, which has no poles in the complex plane, with v replaced by MvM* . 
Then 

V + V* 

M—- — M* 
2 

is also positive definite, and we may conclude that w± vanish on |^| = 1. 

It remains to construct the matrix M. Each of the rows of M separately satisfy 
the same vector Riemann-Hilbert problem as the row vector m, namely each row 
vector of Mvj is regular at Zj. So we need two row vector solutions of this problem, 
with distinct limits at infinity. We start by taking these limits to be given by (1, 0) 
and (0, 1) respectively, so that initially M(oo) = /. These are finite dimensional, 
algebraic problems, so that uniqueness implies existence. But we have already 
proved uniqueness of this problem in Step 2. So there exists a unique M satisfying 
i) and ii) with M{oo) — I. 

We next prove that 

M{z)M* ( -j =M(0). (3.4) 



By the symmetry v- (z) — u*(l/z), we have, in the neighborhood of z 



z. 



■J' 



M{z)M*{l/z) ^ M{z)vj{z)v*{l/z)M*{l/z). 

By assumption the poles of M are such that the row vectors of Mvj are regular at 

z ^ Zj. Now 

and this matrix has the same residue at z = — 1/zj as the matrix w^(z) in (2.18b). 
So v*(l/z)M*{l/z) is also regular in a neighborhood of z = Zj. Since M(z)M*(l/z) 
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has no singularities in the finite z-plane and is bounded at infinity, it is a constant 
matrix, by Louiville's theorem. Letting 2; — > we obtain (3.4). 

Taking z — iin (3.5), we see that M(0) — M{i)M*{i), and so is a positive definite 
matrix. Let A be the positive definite square root of M(0)^^. Replacing M{z) by 
AM{z) we obtain the required gauge transformation satisfying i)-iii) above. D 

4. The Inverse Problem 

In the formulation of the forward problem, the asymptotics of the row vector m 
at z = 0, 00 are given in terms of A(a;) by (2.15). In the inverse problem, A is deter- 
mined by requiring that the solution of the Riemann-Hilbert problem satisfy (2.15). 
In addition, (2.16) must be satisfied, and this leads to an additional constraint on 
A. In the generic case, the constraint (2.17) must also be satisfied, leading to yet 
another constraint on A. One must show that these various determinations of A are 
consistent. 

Consider the matrix Riemann-Hilbert problem 

H^{x,^) =H^{x,^)v{x,^), C £ ^ Hvj regular at z = Zj 

H ^ I as z ^ 00. 

Each row of H satisfies the Riemann-Hilbert problem of the previous section. Hence 
the vanishing lemmas of the previous section imply that this problem has a unique 
solution H when v satisfies the appropriate symmetries. Since detu = dctvj = 1, 
the scalar function dct H has no jumps and no singularities, hence is identically 1. 
Put 



m^ (l,l)Ai7, A- , ^ , 

U A 

Then to, being a linear combination of the rows of H, satisfies the Riemann-Hilbert 
problem, and m ^ (A^^, A) as z ^ 00. We determine A(x) by requiring that (2.15) 
be satisfied, i.e. to. ~ (A, A^^), z ^ 0. (By Theorem 2.6 we may assume that 
v{x,0) = /, hence that H+{x,0) = H^{x,0); hence we may write unambiguously 
H{x,0) for H±{x,0).) This leads to: 

A = X-^Hn{x, 0) + XH2i{x, 0), X'^ = X-^H^ix, 0) + XH22{x, 0); 

hence to two determinations for A: 

2^ gii(x,0) 2^ l-gl2(x,0) 

1-^21(2^,0)' ^ H22{X,0) ■ ^ ■ ^ 

Moreover, (2.16) implies, X^^H^-j^{x,i) + A_ff^(a;,i) = X^^H^2(x, i) + A-ff^(a;,i), 

so 

yL ^ Hiiix,i)- H^2ix,i) ,^ 2) 

H22{x,i)- H+^ix,i)' 
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In the case of reflectionless potentials, or when f = / on the unit circle, we may 
take i7+(x,i) = H^{x,i) in (4.2). In addition, in the generic case, the equations in 
(2.17) imply 

A"^i7+(a;,i) + \H+^{x,i) = 0, \^^H^-^{x,i) + \H:^-^{x,i) = 0; 

hence 



A2 = -±ii^):^ = -lili^. (4.3) 



H+^{x,i) H^i{x,i) 

A similar determination of A is given at —i. 

We must show the equivalence of these various determinations of A. We first 
prove 

Lemma 4.1. Assume the Riemann-Hilbert problem above has a unique solution. 
Then 

RH{x,0)R^H-\x,0), (4.4) 

H+{x,i) = H{x,0)RH_{x,i)R. (4.5) 

Proof. Let H{x,z) = RH{x,—l/z)R. By the s-symmetry of the data, H and H 
satisfy the same Riemann-Hilbert problem. Since the solution of that problem is 
uniquely determined by its behavior at oo, we must have 

H{x, z) = RH{x, -l/z)R == MH{x, z), (4.6) 

where M = M{x) is a matrix independent of z, and detAf(x) = 1. Letting z -^ 
0, oo we find RH{x,oo)R = I = MH{x,0) and RH{x,0)R = M, hence (4.4). 
Writing (4.6) as H{x,z) ~ H{x,0)RH{x,—l/z)R and letting z ^ i from Z^+, we 
obtain (4.5). D 

We are now ready to show the equivalence of (4.1), (4.2), and, in the generic 
case, (4.3). The two expressions in (4.1) are equivalent provided 

Huix, 0) + H2i{x, 0) = 1 - det H{x, 0). 

^From the fact that dctiJ = 1, we see that we must show that Hi2{x,0) + 
H2i(x,Q) = 0; but this follows from the identity (4.4). 

The equivalence of (4.1) and (4.2) reduces to the identity 

Hn{x,0){H^^{x,i) - H+{x,i)) = 

H+{x,t) - H^^ix,i) - H2iix,0)iH+ix,i) - H^^ix,i)). 

(4.7) 

In the case of reflectionless potentials, v{x,^) = /, and H^(x,i) = H^(x,i). In 
(4.7), replace H^{x,i) by H^(x,i), H2i{x,0) by — _ffi2(a;,0); one then obtains the 
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same expression as that obtained by operating on the column vector (1, —1)^ with 
both sides of the matrix equation (4.5). 

In the generic case we write (4.5) in the form H{x, 0)RH+{x, i) = -ff+(x, i)Rv{i) 
and use the expression for v{i) given in Lemma 2.6. We obtain four equations, one 
of which is 

Hii{x,{))H+^{x,i) + Hi2{x,{))Ht2{x,i) = -H+^{x,i). 

This may be used to show the equivalence of (4.1), (4.2). The proofs of the equiva- 
lence of (4.3) with (4.1) in the generic case proceeds along similar lines, and is left 
to the reader. 

Theorem 4.2. Let m{x,z) = {\,l)KH{x,z). Then ■m{x,—l/z) ^ ■m{x,z)R. 

Proof. Since A has been chosen so that m(x,0) — (A, A^^) we have (A, A^^) = 
(l,l)Ai?(a;,0); so 

m{x,z) ={l,l)AH{x,0)H-^{x,0)H{x,z) = {\, X-^)H{x, z) 
^{\\-^)RH{x,-l/z)R 

=(A"\ \)H{x, -l/z)R = (1, l)Ai7(x, -l/z)R 
=m{x,-l/z)R. a 



Theorem 4.3. Let the scattering data v satisfy the symmetry v{^) = Rv ^(^)i?. 
Then |A| == 1. 



Proof. Define H{x, z) — H{x, z). A short calculation shows that 



H+{x,C) = lim H{x,z) = i?-(x,^) 



=H+{x,Ov-HO = H+{x,ORv{OR 
=H-{x,^)Rv{^)R. 

Since the solution of the Riemann-Hilbert problem is uniquely determined by its 
limit as 2; -^ oo and since both H and H tend to / as z — > oo we have H{x, z) 
= RH{x, z)R. Letting z -^ 0+, and recalling that -ff+(x,0) = i7_(x,0), we obtain 
H{x,0) = RH{x,0)R, hence Hii{x,0) = H22{x,0) and Hi2{x,0) ^ H2i{x,0). The 
two expressions for A^ in (4.1) then show that A^ = A^^. D 



For realp, we find that H±{x, z) — H±{x, — z), hence H±{x, i) = H±{x, i). /^From 
(4.2) we sec that A^ is real. We still need to show that this quantity is positive. 
Finally, as in [3] one may prove that H ^ I as x ^ oo, so that (2.14) is satisfied. 
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5. The Initial Value Problems 

We want to solve the initial value problems for (1.2) and (1.3) by the inverse scat- 
tering method. We begin by determining the associated evolution of the scattering 
data. 

The evolution of the scattering data is obtained by standard asymptotic argu- 
ments. /^From [14] the flows (1.2), (1.3) are obtained from the Lax pair 

L = D^-u, u^ikp + q: P = -i— + e(— + k)D 

4 2 

We take e = 2« for the flow (1.2) and e = 2 for the flow (1.3) [14]. The flows are 
obtained from the Lax equations 

L=[P,L] 

when the commutator is evaluated on the submanifold of wave functions of L; that 
is, on functions satisfying Lip + k'^Tp = 0. 
When the Lax equations are satisfled, 

(| - P){D^ -u + E^)4> ={D'-u + i?2)(l -P)4> = 

so 



(|-W 



satisfies (2.1) whenever cf) does. 
Since p ^ to as x ^ ±oo, 



d d 

— P ^ ^ ekD, as a; ^ ±oo. 

dt dt 

so 

^ - P j </- - -efcDe"'-^^ ~ ie kEe'^^'-'^ , as x ^ -oo. 

Since the wave functions are uniquely determined by their asymptotic behavior as 

X -^ — oo, 

(^ - P)^ ^ iekE<i>. 

On the other hand, 0+ — aijj- + bip^, so on S, 

(|-F)0+=(^-P)(a^_+6V+) 

-(a + ikEea)e-'^'' + [b - ikEeb)e'^'' a; ^ oo. 



Therefore, 



d 

(— - P)(/)+ = (a + ikEea)->p^ + {b - ikEeb)i]j+ = iekEiatp^ + bip^) 
dt 
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and 

a^O, b^ 2iekEb. 

Similarly c — —2iekEc. The evolution of the coupling coefficients for the bound 
states, defined by, 

4>{x,t,zf) = cf{t)tl){x,t,zf) 

is derived by similar arguments. We have shown 

Theorem 5.1. Under the flows (1-2) or (1.3) the evolution of the scattering data 
is given by 

r±{x,t,0 = r±(e)e±2^^(^+^'^*); cf{t) = cf{0)e'^<''f\ (5.1) 



where kf and Ej' are the values of k and E at z 



Global existence theorems for each of the flows (1.2), (1.3) are proved by showing 
that the scattering problem can be inverted at all later times, hence, by showing 
that the hypotheses on the scattering data in the two vanishing lemmas are invariant 
under the evolution of the scattering data (5.1). Now, the support of r± and the 
location of the bound states is invariant under the evolution of the scattering data, 
and E is real on S. Therefore, the only remaining constraints to be verifled are 
that efc(^) is real on the support of r±{^) and that iek, E, is real on the bound 
states z^ in the case of real p. 

We restrict the support of r± to lie on the unit circle in the case p real, and on 
the real line for the case p imaginary. Since k is real on the real line, and imaginary 
on the unit circle, and since e = 2i in the Boussinesq case and e = 2 for imaginary 
p, efc(^) is then always real on the support of r±. 

The condition that kfEf be real on the bound states is met for z,- = zw,- since 
in that case, Ej = ismhujj, and kj = icoshcuj. 

The following theorem now follows as a consequence of the invariance of the 
scattering data and the solvability of the inverse scattering problem. 

Theorem 5.2. The initial value problem (1.2) is solvable for — oo < t < oo when, 
at some finite time, the support of the reflection coefficients r±{i) is contained in 
the unit circle, and the bound states are restricted to the imaginary axis and satisfy 
(2.10). 

The initial value problem (1.3) is globally solvable if at some finite time the 
support of the reflection coefficients is contained in the real line, and their are no 
bound states. 
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Data for real p Data for imaginary p 

Global solutions of (1.3) do not exist when there are bound states for the case 
p imaginary. This is illustrated by the example in §6 where we construct a soliton 
for (1.3). In order to prove the vanishing lemma, we need a positivity condition, as 
that of the expression on the right side of (3.5). In the case of real p we used the 
fact that rn{{x^iuj) was real in (3.4). In order to get real wave functions we need 
the coefficients in (1.1) to be real. This means that E is real or imaginary and ikp 
is real, hence that k is real. 

6. Construction of Solitons 

We construct simple solitary waves by solving the Riemann-Hilbert problem in 
the case of reflectionless potentials, that is, when v is the identity matrix and there 
are only discrete scattering data in the problem. When p is real the bound state 
eigenvalues are symmetric with respect to the imaginary axis and the unit circle, 
by Lemma 2.4. In the simplest case we may consider a pair of eigenvalues iuj and 
i/uj, with Lo > 1. 

Our problem is 

(j)j^{x^z) — a{z)ip^{x,z), (f>^{x,iLu) = cip^{x,iLu) (6.1) 

where c is the coupling constant and 

z — iuj 



a(z) = ]— i/'_(x, z) = '0+(x, z ^). (6.2) 

z — i/uj 

We take V'+ == m°^^e'^^^ and assume ?n™* is a meromorphic function in the 
extended z-plane with a single pole at i/cu: 



TO, 



^out 



A-1 + 



w{x) 



z — i/uj 

iFrom (2.4), 

limTO7*(x,z) = A 

z—>0 

hence 



A — A = iujw. (6.3) 
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Invariancc of the wave function under Schwarz reflection in the imaginary axis, 



VIZ. 



leads to the constraint 



V'+(X, -Z) =-0+(.T,z), 



z — IJUJ z — i/uj 



for aU z. Hence we must have 



A = A, w — —w. 



Setting 



W=-—-, A = e^/2^ Lo = e", P ^ -iEiiu) = sinh a, 

zip 



we obtain from equations (6.1-6.3) the relations 



P 

W — — a;/3sinh— . 
2 



Solving for P, we obtain 



e^ = l 



2c 



/3w3 exp (2/3a;) + cw/3 



Then, from (2.5), 



'4-;l(l-'-'"^ 



4 dx 



P 



2/3c 



2 ^2 exp (2/3a;) + c 



(6.4) 



(6.5) 



By Theorem 5.1 the time evolution of the coupling constant c is 

c(i)=coe2*^'"h2a^ 

The eigenvalues for imaginary p appear in fours, with the symmetries given in 
Theorem 2.3. We construct a soliton corresponding to four bound states with purely 
imaginary z^ ~ ilo, hence 



Then 



iw, i/w, —i/oj, —iui. 
z — icu z + i/uj 



a{z) = 



z + iuj z — i/uj^ 
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and 

(f)-f-{x, z) — a{z)ip^{x, z), ip^{x, z) = xp+{x,z). (6-6) 

^From Theorem 2.3 

(f)^{x,iLo) = cil)^{x,iLo), (j>^{x, —i/uj) ^ cxp+{x,—i/Lo) (6-7) 

where c is the couphng coefficient. 

By Theorem 5.1, the time evolution of c is given by 

c{t) = c(0)e2^^'^-^* ^ c(0)e-2^(^'"'^ 2a)t^ 

Now ip+{x, z) = TO™* (a;, z)e*^^. We look for a wave function of the form 

TOr(x,z)=A-i+^^^+^^. 

The symmetry 

to7*(x, z) ^ m^"*(x,-z-i) 

implies that 

out/ ^ T-1 2Wi ZW2 

m_f (x, z) = A 



1 + icijz 1 — iz/u 

The requirement that to,^"*(x,0) — \{x) implies that A = A^^, hence that |A| = 1. 
We then have 

A-A-=f^^ + -^Vf^V + T^^V (6-8) 

\z + luj l — iz/ujj \Z — l/UJ 1 + lUJZj 

Since the right side is independent of z, 

w\ + ijli^W2 — 0; 
and 

TOr(.,.)^A- + ^ + -fli:^. 

z + luj lujyi + lujz) 
Equations (6.6) and (6.7) now lead to an equation for w ~ wi: 

w = « (a-^ + -^ - ,. ,^ , ) e-^^^^"^ " (6.9) 

\ ZiLO 2iL0^ smh a/ 

where 

e" = Lu, K = 2iajctanh a. 

Taking complex conjugates of (6.9) we obtain the following system of equations 

for w and w: 

a b\ f w\ _ fa 
b a I \ w j \ a 
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where 

a=l + -sechae-2^ '''"'' ", b ^ -ctanh ag-^^'^'"^ ", 

CO 

and 

(J = -2icjcAtanh ae'^'"''''^ ". 

The solution of this 2x2 system is 

acr-ba 
w = — - — , (6.10) 

where A = |ap — |6p. To determine whether a global solution exists, we must 
evaluate A with c(i) given as above. 
Taking z = in (5.8) we find 

3A=^^ = -^^^. (6.11) 

Equations (6.10), (6.11) arc coupled equations for A and w. 

We note that the solitons for imaginary p always have poles on the real axis. 
These occur at the zeros of A, hence, whenever \a\ = \b\. To show that A always 
vanishes, first set C =■ coe^^'^^'"'^ ", so that < |C| < oo. As t varies, |6| — 
jCtanh a\ is constant, while \a\ ranges over the interval (1 — |C|e^"secha, 1 + 
|C|e~"sechQ!). Therefore A vanishes whenever (we may assume a > 0) 

1 - |C|e^"sechQ! < |C| tanh a < 1 + |C|e""secha; 

this inequality simplifies to 

0<(e2" + l)(l-l/|C|) <4, 

which is always satisfied for a range of \C\ in (0, oo). Hence the soliton always has 
poles for real x and t. 

7. Conclusions 

We have proved the existence of global solutions of the initial value problems of 
the flows (1.2) and (1.3) under complementary conditions on the scattering data. 
Two questions are often raised about these results. Firstly, can the conditions on the 
initial data be expressed in terms of the initial values u(x^ 0) and w(x, 0)? Secondly, 
are the conditions sharp or merely sufficient conditions for global existence? 

As to the first question, it is typical in the theory of the scattering transform 
that the picture is often simpler on the scattering side than on the potential, or 
spatial side. This is perfectly analogous to the features of the Fourier transform. 
For example, under the Fourier transform, convolution, a non-local operation, is 
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mapped into pointwisc multiplication, a local operation. In the case of completely 
integrable systems, action-angle variables are expressed very conveniently on the 
scattering side. [NMPZ]. Thus, without saying that simple conditions for global 
existence cannot be placed directly on u and w, we should also not be surprised 
that the solvability conditions are expressed so simply in terms of the scattering 
data. 

The second question is in fact one we wish to address in future work, but we 
conjecture that the conditions are in fact sharp. One already sees in the example 
of the soliton in §6 for the case of imaginary p: it has poles for all time. We 
conjecture further that if one takes as initial data reflection coefficients r±{^) with 
support on the wrong section of S, then singularities will develop in finite time. For 
example, we conjecture that if for realp one chooses for r-i-(^) smooth functions with 
compact support on the real line, singularities will develop in finite time. It would 
be interesting to use the Riemann-Hilbert formulation to investigate the nature of 
the blow up of the solution. Presumably it would develop a pole, since that is the 
nature of the singularities of the solution of the inverse scattering problem. 

The situation is in some respects similar to that in the AKNS system, in which 
some kind of symmetry on the potential is needed in order to guarantee global solv- 
ability. For example, the nonlinear Schrodinger equation is obtained for Hermitian 
or skew-Hermitian potentials, the modified KdV equation is obtained for real sym- 
metric or skew symmetric potentials. The symmetry of the potentials is related 
to the positivity of the quadratic conservation laws. Without the symmetry, these 
conservation laws are not positive definite; and global solutions of the unsymmetric 
AKNS systems do not in general exist (cf. the discussion in [2].) 

Appendix 

Proposition Al. Let p, q G i"'^(K). Then for z E lAi the wave function -0+ can 
he written in the form 

^+{x,z)^(Aa{x,E) + -^Ai{x,E)\e'^'\ zeUi, (Al) 

where Aq, Ai arc analytic in "^E > 0, continuous onto the real E-axis, and 

Ao-A"\ Ai-0, £;^oo,3£;>0. 

The form of the wave functions in the other domains can be obtained from the 
s symmetry. 

Proof. Using the identities 

^1 I 2E , „ 1 

z = 2E--, — = 1, k = E--, 

z 
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and the ansatz (Al), the integral equation (2.3a) can be written as the eoupled 
system of Volterra integral equations 

Ao{x,z)=l- / — {{iEp{y) + q{y))Ao{y,z) + ipAi{y,z)) dy, 

Ai{x, z) = - / — — iiqiy) - 2iEp{y))Ai{y, z) - ipAo{y, z)) dy. 

Jx 2ii; 

These may be solved by successive approximations for "^E > in the standard 
way, and the analyticity of Aq and Ai established. Letting _B ^ oo we obtain the 
limiting equations for Aq and Ai, 

Ao{x) = l-- p{y)Ao{y)dy, Ai(x) = p{y)Ai{y)dy, 

from which the result follows. D 

Lemma A2. The wave function ^p^{x,z) has the representation 

ij+{x,z) = A U'^- + £ {K^{x,y) + -^K,{x,y))e''y)d^ . (A2) 

Similar representations hold for the other wave functions tp^ , (f>± . Here, Kq and Ki 
are smooth kernel functions, with properties similar to the kernels in the Gel fand- 
Levitan-Marchenko integral equations in the scattering theory of the Schrodinger 
equation. 

Remark: Kaup, ([10], (3.7)) assumes the representation (A2) in his derivation of 
the GLM equations. 

Proof. These Fourier representations are proved by arguments similar to those used 
in the standard GLM theory. Since Aq and Ai are analytic in the upper half E- 
plane, and since Aq — \^^ -^ a.s E ^ oo, they can be represented, for each x, as 
Fourier transforms 

/o /-O 

Hoix,s)e-'^''ds, Ai= Hi{x,s)e-'^''ds, 

Writing 

r" 

e*-^^Ao = Ae*-^^ + / Ho{x, s)e'^^''~'Us, 

J —OQ 

putting y = X — s^ and changing variables in the integral, we obtain 

Aoe*^- - A ("e*^- + y" Ko{x, y)e'^ydy] , 

where KQ(x,y) = X^^Hq(x,x — y). The other term is dealt with in the same 
way. D 
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